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Gbrtler  instability  for  boundary- layer  flows  over  generally  curved 
walls  is  considered.  The  full  linearized  disturbance  equations  are 
obtained  in  an  orthogonal  curvilinear  coordinate  system.  A  perturbation 
procedure  to  account  for  second-order  effects  is  used  to  determine  the 
effects  of  the  displacement  thickness  and  the  variation  of  the  stream¬ 
line  curvature  on  the  neutral  stability  of  the  Blasius  flow.  The  pres¬ 
sure  gradient  in  the  mean  flow  is  accounted  for  by  solving  the  non¬ 
similar  boundary- layer  equations.  Growth  rates  are  obtained  for  the 
actual  mean  flow  and  compared  with  those  for  the  Blasius  flow  and  the 
Falkner-Skan  flows.  The  results  demonstrate  the  strong  Influence  of  the 
pressure  gradient  and  the  nonsimilarity  of  the  basic  flow  on  the 
stability  characteristics. 
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I.  Introduction 

Centrifugal  instability  in  boundary-layer  flows  over  concavely 

curved  walls  was  first  demonstrated  theoretically  by  Gortler^.  The 

disturbances  are  in  the  form  of  counter-rotating  streamwise  vortices, 

2 

called  Gortler  vortices.  Wortman  gave  a  detailed  flow  visualization  of 
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these  vortices  on  slightly  curved  walls.  Bippes  conducted  experiments 
on  walls  with  radii  of  0.5  and  1.0m  so  that  the  generated  vortices  were 
fairly  strong. 

These  vortices  by  themselves  may  not  lead  to  the  transition  of 
laminar  flows  to  turbulent  flows.  The  influence  of  steady  streamwise 

vortices  on  two-dimensional  Tol Imien-Schl ichting  waves  was  studied 
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experimentally  by  Aihara  ,  Tani  and  Sakagami  ,  and  Tani  and  Aihara  . 

They  concluded  that  the  Gortler  vortices  indirectly  affect  the  transi¬ 
tion  by  inducing  a  spanwise  variation  in  boundary-layer  thickness,  at 
least  when  the  radii  of  curvature  are  not  extremely  small.  After  re¬ 
viewing  the  experimental  work  on  the  subject,  Nayfeh7  presented  a 
theoretical  model  for  studying  the  effect  of  streamwise  vortices  on 
Tol Imien-Schl ichting  waves.  He  showed  that  such  vortices  have  a  strong 
tendency  to  amplify  three-dimensional  Tol Imien-Schl ichting  waves  having 
a  spanwise  wavelength  that  is  twice  the  wavelength  of  the  vortices. 
Therefore,  the  study  of  Gortler  instability  is  of  importance  not  only  in 
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heat  exchange  problems  but  also  In  the  analysis  of  the  transition 

region  over  concavely  curved  walls. 

Improvements  and  extensions  of  the  Gbrtler  model  were  considered  by 
8  Q 

many  authors.  HSmmerlln  relaxed  some  of  the  assumptions  made  by 

g 

Gortler  regarding  the  treatment  of  the  curvature  terms.  He  accounted 

for  the  variation  of  the  streamline  curvature  normal  to  the  wall  by 

using  a  flow  that  has  the  properties  of  boundary- layer  flows  along  a 

wavy  wall  in  the  valley  positions.  For  these  flows,  the  streamline 

curvature  decays  exponentially  away  from  the  wall.  Tobak^  considered 

the  case  of  a  wall  having  a  finite  extent  of  streamwise  curvature.  The 

effect  of  this  curvature  in  the  final  equations  appears  as  a  multiplier 

to  the  Gortler  number.  Herbert^1  gave  a  comprehensive  survey  of  the 

effort  of  these  authors  and  obtained  very  accurate  numerical  solutions 

to  their  models. 

1  ? 

Smith  Incorporated  the  transverse  component  of  the  basic  flow  and 

some  of  the  higher-order  curvature  terms.  Furthermore,  he  considered 

spatially  rather  than  temporally  growing  vortices.  Smith's  viewpoint  is 

consistent  with  the  experimental  observations.  Smith  used  body-oriented 

axes  and  a  Taylor-series  expansion  of  the  metric  coefficient  for  small 

distances  normal  to  the  wall.  He  kept  two  terms  In  the  expansion.  A 

complete  survey  and  critical  discussion  of  the  role  of  the  coordinate 
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system  in  the  Gortler  problem  was  given  by  Floryan  and  Saric  .  They 
reformulated  the  problem  In  an  orthogonal  curvilinear  coordinate  system 
constructed  from  the  potentlal-and  stream-lines  of  the  Invlscld  flow 
over  the  curved  wall.  They  pointed  out  how  this  coordinate  system 
overcomes  the  singularity  In  the  metric  coefficient  associated  with 
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body-oriented  axes.  They  considered  the  case  of  a  circular  arc  and 
obtained  the  basic  approximations  of  the  mean-flow  and  disturbance 
equations.  They  concluded  that  the  correct  zeroth-order  base  flow  for 
the  stability  analysis  is  the  Blasius  flow.  They  based  the  scaling  of 
the  disturbance  in  the  transverse  component  on  a  viscous  scale,  namely 
v/6r,  where  v  is  the  kinematic  viscosity  and  6r  is  a  reference  boundary- 
layer  thickness.  As  a  result  of  this  ordering,  the  transverse  component 
of  the  basic  flow  has  a  first-order  effect  on  the  stability  analysis. 
This  fact  was  also  pointed  out  by  Herbert11  after  introducing  the 
Gortler  number  into  Smith's  equations.  Moreover,  the  streamwise  varia¬ 
tion  of  the  transverse  component  of  the  basic  flow  appears  also  in  the 
first-order  equations.  This  term,  which  has  a  strong  Influence  on  the 

stability  characteristics  (in  particular  the  neutral  stability  curve), 
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is  missing  from  the  analyses  preceding  that  of  Floryan  and  Sarlc  . 

Here,  we  are  interested  in  flow  over  walls  with  general  distri¬ 
butions  of  curvature.  Specifically,  we  consider  the  effect  of  the 

streamwise  pressure  gradient  on  the  stability  characteristics.  Follow- 
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ing  Floryan  and  Saric  ,  we  employ  a  coordinate  system  based  on  the 
potential-  and  stream-lines.  The  inviscid  flow  and  the  metric  coeffi¬ 
cients  of  the  coordinate  system  are  obtained  by  the  method  of  Davis1^ 
for  constructing  body-fitted  coordinate  systems.  His  method  is  based  on 
the  Schwarz-Christoffel  transformation  generalized  for  curved  surfaces. 

In  Section  II,  the  two-dimensional  basic-flow  equations  are  derived  in  a 
general  orthogonal  curvilinear  coordinate  system.  The  obtained  equations 
are  specialized  for  body-oriented  axes  and  for  potential-  and  stream¬ 
line  coordinates.  In  Section  III,  the  complete  linearlized  Navier-Stokes 


equations  are  derived  for  the  disturbance  quantities.  The  quasi-parallel- 
flow  assumption  and  the  normal-mode  solution  are  discussed.  Then  a 
method  of  solution  is  presented.  Terms  of  order  smaller  than  0(1)  are 
treated  as  perturbation  quantities  to  a  basic  system  that  includes  the 
0(1)  terms.  The  method  of  solution  of  the  basic  system  is  presented  in 
Section  IV.  A  new  method  of  applying  the  asymptotic  boundary  condition 
far  from  the  wall  is  presented.  In  Section  V,  a  perturbation  procedure 
is  presented  to  account  for  the  small  quantities  in  the  disturbance 
equations.  Results  and  their  discussion  are  presented  in  Section  VI. 

II.  The  Basic  Flow 


We  consider  steady  two-dimensional  viscous  incompressible  flows 

over  a  curved  wall.  The  fluid  properties  p*  and  p*  are  assumed  to  be 
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constants.  Following  Floryan  and  Saric  ,  we  employ  an  orthogonal 
curvilinear  coordinate  system  <J>*,  ^*,  and  z*.  Later,  we  identify  4>*  and 
4>*  as  the  potential  and  stream  functions  of  the  two-dimensional  inviscid 
flow  over  the  given  wall,  respectively,  and  z*  as  a  rectilinear  coordi¬ 
nate  normal  to  the  plane  of  the  basic  motion,  see  Fig.  1.  The  Navier- 
Stokes  equations  in  this  coordinate  system  are  given  in  Appendix  A.  We 
introduce  the  following  nondimensional  variables: 


=  <J>*/U*L*,  ip  =  ip*/U*L* 


h.  =  h*U*,  h.  =  h*U* 

<p  <j>  oo  Ip  ip  co 


=  u*/U*,  v  =  v*/U*,  w  =  w*/U*, 


p  =  p*/p*U*2 


We  define  a  characteristic  Reynolds  number  Re  =  P*U*L*/v*»  where 
p*,  and  p*  are  the  free-stream  values  of  velocity,  density,  and 


(1) 

(2) 

(3) 

u*. 
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viscosity  and  L  is  a  characteristic  length,  the  distance  from  the 
leading  edge  of  the  curved  wall  to  the  position  of  maximum  height,  see 
Fig.  2. 

The  equations  of  Appendix  A  for  the  two-dimensional  case  (l.e., 
w  =  0,  3/3z  =  0)  become 


^-Momentum  Equation 


U  3u  V__  3U  _  V2 

F7  34>  h,  3ip  h,h, 
<f>  V  <fi  ^ 


.  uv  .  1_  9£  _  1_  _J _ f~3  / ?  3u 

3<j>  h^  dip  h0  3<j>  ~  Re  h^  jj<p  ^  d<p 


+  2  S  ^  *  h [lv  k  +  (V]  +  b  (V 


h,  ^  ,,  3h,  ,  3h  3h,~) 

+  _£  <L  (iL.n  _JL  _  21—  —  +  — _ Jh  [ 

\  3*  vJ  3*  \  3*  3*  3*  J 


(4) 


ip-Momentum  Equation 


u_  v_  3v  uv  w"tp  _  u ‘  w"<t>  1_  3£  1_  J _ I  jL  r jfc  3_  ,v_^ 

K  ^  hw,  ^  3*  ’  hA,  W  h,i,  W  ~  Re  hA,  W  Lh*  94>  A, 


3h,  .,2  3h 


+  If  *  k  <2  ^  U  *  2  5*  si*’  +  [hjl*  (V 

*  \  I*  lV]  ** ' 2[*A  *  «] 


*P 


O 

3_ 

3<j> 


(5) 


Continuity  Equation 


%  <y) +  <y> =  0 


(6) 


We  assume  that  the  metric  coefficients  h.  and  h.  and  their  derivatives 

<p  v 

that  appear  in  Eqs.  (4)-(6)  are  all  0(1)  or  less.  Hence,  regions  of 
severe  or  discontinuous  curvature  are  excluded. 

Following  Prandtl,  we  introduce  the  stretched  variable 
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_  1  / 
ip  =  Re  ' 2 


(7) 


i  / 

Then,  It  follows  from  Eq.  (6),  that  v  must  be  0(Re  /2)  In  the  least- 
degenerate  limit.  Hence,  we  let 


v  = 


v 


(8) 


where  v  =  0(1).  Substituting  Eqs.  (7)  and  (8)  Into  Eq.  (4)  and  neglect¬ 
ing  terms  0{Re_1)  we  obtain,  in  terms  of  the  unstretched  variables. 


u  Su  ,  v  3u  ,  uv  ,  1  9p  1  1  ( 9_  r ^  9_ 

%  \  **  Vv  \ " Re  ^vi3*  \  ^ 


h,  a  it  3hM  1 

+  JiL.  /«_)  _i  C 
h*  *  \  *  J 


(9) 


Similarly,  the  ^-momentum  equation  reduces  to 


u2 
TT 
<P 


3i(i 


l£ 

dip 


The  continuity  equation  is  unchanged. 


(10) 


A.  Body-Oriented  Axes 

For  a  body-oriented  axes  we  Identify  <j>  as  x,  the  arc  length  mea¬ 
sured  along  the  body  contour,  and  ip  as  y,  the  normal  to  that  contour. 
Hence 


\  =  1  +  k(x)y 


(ID 


(12) 


where  k(x)  is  the  local  curvature  of  the  contour.  Now,  Eqs.  (9),  (10), 
and  (6)  become 


TTJy 


kuv 

l~TTy  +  r 


__  3£ 
ky  3x 


!_  r32U  +  k  9u_ 

Re  LSyT  f  +  Icy  9 y 


ku2  _  5p 
1  +  ky  "  3y 


(14) 


|i  t  !„  [(,  +  ky)V]  =  0  (15) 

liquations  (13)-(14)  have  been  derived  by  Schultz-Grunow  and  Breuer^, 
among  many  others  (see  Van  Dyke^  for  a  complete  review),  as  an  appro¬ 
priate  set  of  equations  to  describe  boundary-layer  flows  with  second- 
order  curvature  effects. 


B.  Equations  in  Potential  and  Stream-Lines 

Here,  we  identify  <p  and  ip  as  the  potential  and  stream  functions  for 
the  inviscid  irrotational  flow  over  the  given  wall.  Furthermore,  we 
choose  \p  =  0  to  coincide  with  the  body  surface  and  $  =  0  to  coincide 
with  the  leading  edge.  It  follows  from  the  definitions  of  h^  and  h^ 
that 


hj,  =  h^  =  1/U(M)  =  h  (16) 

where  U  is  the  magnitude  of  the  inviscid  velocity.  Actually,  it  is  the 
4>- component  of  the  velocity  since  the  velocity  normal  to  a  streamline  is 
by  definition  equal  to  zero.  From  the  irrotational ity  condition  we 
obtain 


1  3h 


and  from  the  incompressibility  condition  we  obtain 


k  =1^ 
<p  h7  d<p 


(17) 


(18) 


where  k,,  and  k,  are  the  local  values  of  the  curvature  of  a  streamline 
and  an  equipotential  line,  respectively. 
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P  J  in  ■  wiW 


In  this  coordinate  system,  Eqs.  (9),  (10),  and  (6)  become 


.  3u  ,  .  3u  .  3h  ,  3d  .  1  32u 

hu  3<J>  hv  3ip  UV  3t|;  "  _h  if  Re  itp" 


(19) 


u2  3h  _  3p 
h  dip  _  a* 


(20) 


k  (hu)  +  If  (hv)  -  0 


(21) 


where  the  term  ^  was  neglected  from  Eq.  (19)  since  It  is  0(Re“1). 
consistent  with  our  assumption  that  the  metric  coefficients  and  their 
derivatives  are  <  0(1).  Equations  ( 1 9)-( 21 )  are  alternatives  to  Eqs. 
(13)-(15)  of  the  body-oriented  axes. 

We  identify  two  types  of  curvature  effects  on  the  basic  state.  The 
first  is  a  viscous  effect  that  manifests  Itself  through  the  terms 
uv9h/3i|;  In  Eq.  (19)  and  (u2/h)3h/3ip  in  Eq.  (20).  As  shown  by  Van  Dyke, 
this  effect  is  a  local  one.  The  second  is  an  Inviscid  effect  that 
manifests  itself  in  the  pressure  gradient  term  3p/3<|>.  Since  in  subsonic 
flows  this  effect  is  a  global  one,  it  depends  on  the  distribution  of  the 
surface  curvature  everywhere.  Consider  for  example  the  flow  over  the 
hump  shown  in  Fig.  2.  The  central  part  of  the  hump  has  an 
appreciable  curvature,  specifically,  a  curvature  of  0(1),  while  the 

_i  / 

upstream  part  of  the  hump  has  a  negligible  curvature,  say  0(Re  '2). 
Nevertheless,  the  pressure  gradient  3p/3(J>  In  the  latter  region  may  not 
be  negligible.  One  of  the  main  purposes  of  the  present  paper  is  to 
evaluate  the  influence  of  this  pressure  gradient  on  the  stability 
characteristics  of  Gortler  vortices.  Hence,  for  this  purpose,  we 
neglect  the  viscous  effect  of  curvature  on  the  basic  state  because  it  is 
a  local  effect. 
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We  introduce  the  new  variables 


u  =  hu 

(22) 

\l  =  hv 

(23) 

Substituting  Eq.  (22)  and  (23)  into  Eqs.  (1 9)- (21 )  and  neglecting  terms 

,-V  2 


of  order  Re 


u  M  +  v  M 

u  3<J>  3i|j 


we  obtain 

(u 


1  3h  /~2 


h  34> 


i)  +i_  afo 

1 '  Re  34^ 


(24) 


(25) 


3u  +  3v 
34>  3t|> 


=  0 


(26) 


The  relation  h2  =  -r  which  is  valid  for  the  inviscid  flow,  has 

d4>  n  34> 

been  used. 

Now,  we  introduce  the  Levey-Lees  variables 


£  =  4>  (27) 

n  =  ir(Be  (28) 

/2£ 

F  =  u  (29) 

V  =  Vi^Re  /2£  -  nF  (30) 

Then  the  4>-momentum  and  the  continuity  equations  take  the  standard  form 
2£FF^  +  VFn  +  Bo(F2  -  1 )  -  F^  =  0  (31) 

2£Fr  +  V  +  F  =  0  (32) 

£  n 

where 


Bo 


3h 
h  34> 


(33) 
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where  H*  is  the  maximum  height  and  B  is  a  constant.  The  curvature 
variation  of  this  wall  is  shown  in  Fig.  3  as  a  function  of  x*/L*. 
Equation  (38)  shows  that  y*  decays  exponentially  as  x*  -*■  ±°°.  For 
moderate  values  of  B  (about  10),  y*/H*  ~  2e"10  at  x*  =  0,  which  is  very 
small.  Therefore,  for  the  purpose  of  calculating  the  Invlscld  flow  over 
the  hump,  we  assume  that  the  wall  Is  specified  by  Eq.  (38)  In  the 
interval  0  <  x*/L*  <  2  and  it  coincides  with  the  x*-axis  outside  this 
interval . 
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Davis  developed  an  accurate  numerical  method  for  integrating 

Schwartz-Christof Tel  transformations  generalized  for  curved  surfaces. 

We  employ  Davis'  method  to  find  a  function  f(c)  which  maps  the  upper 

half  of  the  c-plane  onto  the  region  above  the  curved  wall  in  the  z-plane 

A  uniform  flow  parallel  to  the  £-axis  in  the  s-plane  is  transformed  into 

a  flow  over  the  hump  in  the  physical  plane.  Far  from  the  hump  the  flow 

14 

returns  to  a  uniform  stream.  The  paper  by  Davis  contains  sufficient 
information  for  programming  the  method.  However,  application  of  Davis' 
method  to  the  wall  treated  in  the  present  work  is  included  In  Appendix  B 

D.  A  Method  For  Solving  The  Basic  Flow  Equations 

Once  Bo  is  determined  from  the  inviscid  flow  solution,  then  the 
solution  of  Eqs.  (31)  and  (32)  subject  to  the  boundary  and  initial 
conditions,  Eqs.  (35)-(37),  is  obtained  by  using  a  finite-difference 
marching  technique.  A  second-order  accurate  Scheme  is  employed.  A 
Newton-Raphson  process  is  used  to  quasi-1 inearize  the  nonlinear  terms  in 
the  momentum  equation.  Hence,  the  momentum  and  continuity  equations  are 
coupled  in  their  linearized  forms.  This  is  known  as  the  Davis  coupled 
scheme  (see  Blottner^7).  The  ^-derivative  is  replaced  by  a  3-point 
backward  formula  while  the  n-derivatives  are  represented  by  central 
differences.  At  the  first  two  streamwise  locations,  the  profiles  are 
taken  to  be  the  Blasius  profile. 

For  the  Falkner-Skan  flows,  B0  is  a  constant  and  a  self-similar 
solution  for  Eqs.  (31)  and  (32)  can  be  obtained.  Since  the  equations 
are  ordinary-  rather  than  partial-differential  equations,  we  use  a 
fourth-order  Runge-Kutta  method  coupled  with  a  shooting  technique  to 
satisfy  the  two-point  boundary  conditions  (35)  and  (36).  The 
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asymptotic  condition  (36)  is  satisfied  at  n  ■  18,  which  proved  to  be 
sufficiently  large.  The  step  size  An  Is  taken  to  be  0.01  for  all  the 
numerical  results  presented  in  this  paper. 


The  basic-flow  quantities  that  appear  as  coefficients  in  the 


perturbation  equations  of  the  next  section  are 

hu  =  F 

(39a) 

9hu  _  9F  .  9n  _  y'Re  9F 

9i p  9n  9ty  9n 

(39b) 

hv  =  *-±JlF 

(39c) 

/2Re<j> 

m  =  Ihv  .  9ii  a  1  9_  (v  +  F) 
dip  9n  9^  2<j>  an  [  ; 

(39d) 

=  Ihv  9£  +  9hv  M  =  _J _ 1  r.  (V  +  nF\ 

9<()  95  9<j>  9n  3d>  2v^Ie  $  2  L 

-  n  (V  +  nF) 

+  25  ~  (V  +  nF)] 

(39e) 

For  further  reference  we  define 

U0  =  hu  =  F 

(40a) 

p  -  9hu  _  1  9  F 

0  "  /r¥  **  "  /2  3ri 

(40b) 

V0  =  hv  =  —  rlF 

J2 

(40c) 

M  _  .  9hv  _  1  9Vo 

Ho  =  ♦  dip  '  sz  9n 

(40d) 

r  -  /Da  !  2  -  1  U  ^Hp  ,  £.  9Vq 

G°  =  *  9* - I  V°  '  +  «  95 

(40e) 

The  behavior  of  the  basic-flow  solution  as  i|>  »  is  also  required.  It 


follows  from  Eq.  (36)  that 

U0  -  1 

(41a) 

E0  -  0 

(41b) 

Equation  (32)  can  be  rewritten  as 

IK.  +  i  =  (■)  _  p)  +  25  l  ~ 
9q  K  1  4  95 


which  upon  integration  with  respect  to  n  yields 

V  +  n  =  o/n  0  -  p)dn  +  25  ~  P  (1  -  F)dn 

The  integral  in  this  equation  is  finite.  As  n  -*•  “  it  Is  proportional  to 
the  displacement  thickness,  hence 

V  +  n->/2g(5)  as  n  “ 

where  g(5)  is  a  function  of  5.  For  the  Falkner-Skan  profiles,  g(5) 
becomes  an  absolute  constant  that  depends  on  Bo-  For  the  Blasius  flow 


constant  is  0 

.8604. 

Therefore, 

Vo  ->  g(5) 

as 

Ip  -*■  CO 

(41c) 

Ho  -*  0 

as 

\p  -*  CO 

( 4 1  d ) 

Go  g(£) 

+  5  ^ 
5  d5 

(41e) 

It  is  worth  noting  that  the  transverse  component  of  velocity  V0  and 
its  streamwise  derivative  Go  do  not  vanish  outside  the  boundary  layer. 
This  is  because  the  displacement  thickness  effect  is  neglected  in  the 
theory  thus  presented.  However,  we  consider  the  flow  due  to  the  dis¬ 
placement  thickness  and  its  influence  on  the  stability  characteristics 
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only  for  the  Blasius  flow.  Also,  we  note  that  the  transverse  component 
of  the  boundary-layer  solution  in  the  present  coordinate  system  does  not 
grow  indefinitely  as  »  for  a  non-zero  pressure  gradient.  This  is  In 
contrast  with  the  behavior  of  the  v-component  In  the  body-oriented  axes 
where  it  grows  linearly  with  y  as  y  -*•  ».  Therefore,  we  Identify  the 
present  coordinate  system  as  optimal  to  first  order.  Had  we  used  the 
stream  lines  and  potential  lines  of  the  Invlscid  flow  over  the  displace¬ 
ment  body  we  would  have  obtained  optimal  coordinates  to  second 
18  19  20 

order  *  ’  .  In  such  coordinates  the  behavior  of  the  boundary- layer 

solution  produces  the  inviscid  outer  flow  to  second  order  (l.e..  It 
includes  corrections  due  to  displacement  thickness). 

E.  Effect  of  the  Displacement  Thickness  on  the  Blasius  Flow 
1 8 

Van  Dyke  showed  that  the  outer  solution  of  the  stream  function 
I*0  for  the  second-order  Blasius  flow  is 


4'°  =  y  -  Real  (A  +  iy)  (42) 

/Re 

and  the  inner  solution  is 

T1  .  ^  .0  (43) 

/Re  m  Re 

where  x  and  y  are  Cartesian  coordinates  (they  are  Identical  to  <p  and  ip 
for  the  case  of  a  flat  plate),  fi(n)  is  the  Blasius  profile  and  Bi  = 
1.21678.  A  second-order  composite  expansion  can  be  formed  from  and 
vp1  according  to 

r  i  0  4  0 

Y  =  y  +  4*  -  (4*  ) 

or 


15 


r  =  —  <fi  +  6i[l  -  ReaHl  +  1  *■  )  /2j  r  (4* 

/Re  L  J 

Now,  the  v-component  and  its  x-derivative  defined  by  Eqs.  (40c)  and 

(40e)  become 

V02  •  V01  -  ^-  +  Real (1  +  i  J)"V2  (4£ 

Z2  /2  x 


Go 2  =  Goi  +  Real  [— — -1  — 

2/2  2/2  (1  +  iy/x) 


l  '2 


where  V02  and  602  include  displacement  effects  while  V0i  and  G0i  do 
not.  We  note  that 

Vqi  8i  1/2  and  Goi  -*■  -8i/2  Z2 


as  y  ■>  «,  and  hence  V02  and  G02  vanish  as  y  +  as  they  should. 


III.  The  Disturbance  Equations 

We  consider  the  stability  of  a  basic  flow  over  a  two-dimensional 
curved  wall  with  respect  to  three-dimensional  disturbances  in  the  form 
of  alternating  streamwise  vortices,  known  as  the  Gortler  vortices. 
However,  we  need  not  specify  this  type  of  disturbance  at  this  stage  of 


analysis.  Instead  we  assume  the  total  flow  in  the  form 
u  =  u0(<M)  +  ui(<J>,^,z) 


v  =  v0(<p,V)  +  vx  (4>,vp,z) 


wi(<M,z) 


P  =  Po(<M)  +  Pi(<M,z) 


(47a) 


(47b) 


(47c) 


(47d) 


where  <p  and  t  are  the  potential  and  stream  functions  of  the  Invlscld 
basic  flow  and  z  is  a  rectilinear  coordinate  normal  to  the  plane  of  the 
basic  flow;  u,  v,  and  w  are  the  velocity  components  In  the  t»  and  z 
directions,  respectively;  and  p  Is  the  static  pressure.  The  subscript 
'O'  denotes  a  basic-flow  quantity  while  the  subscript  T  denotes  a 
perturbation  quantity.  Substituting  these  total  flow  quantities  Into 
the  Navier-Stokes  equations  (see  Appendix  A),  subtracting  the  basic-flow 
quantities,  and  keeping  linear  terms  in  the  perturbation  quantities,  we 
obtain 


t-Momentum  Equation 


~  3u,  ~  3ui  .  r  3un  .  Uo  ah  Vn  ah  1  V2h-, 

Vo^"Uoat  +  ["at  h  at  -  ^  at 


,du0  2v0  3hv  2^1  3h  dv±  ,  2_  1  ah  av^  .  3 

'■3\p  h  3t;  1  "  Re  h  a*  at  Re  h  ty  3*  "  at 


t-Momentum  Equation 


Re  V*Vl  ”  Vo  1^ 
1  V2h 


3vi 

at 


-  u0  -Kr1-  +  [- 


3vo 

at 


uc  ah  vt.  ah 

'  h  at  h  at 


1 3v  r 


■_  yyn  / _ 

'  Re  h  JVl  "  vat  ‘  h  at 

2_  i  aji  au^  ,  a^  _  n 

"Re  h  at  at  a,lj  ” 


2u^  ah\  ,  2__  i_  ah  au i 
;Ul  Re  h  at  at 


at 


at 

(48) 


(49) 


z-Momentum  Equation 


1 

Re 


V2Wi  -  v0 


awi 

at 


Continuity  Equation 
ahv 


ahui 

at 


at 


+  h 


2  9Wl  =  0 

az  u 


(50) 


(51) 


where 


u0  a  hu0,  Vo 

=  hv0. 

(52) 

V*  =  14  +  14 

v  9^  3i Y 

+  h2 
h  az7’ 

(53) 

*  ★  ★  ★ 

(54) 

Re  -  p  UJ.  /u 

Equations  (48)  -  (51)  are  the  complete  linearized  form  of  the  Navier- 
Stokes  equations. 

A.  Scaled  Variables 

The  disturbance  equations.  Eqs.  (48)  -  (51),  are  non-dimensional ized 

★ 

with  respect  to  an  arbitrary  length  L  .  However,  for  stability  analysis 
of  boundary  layers,  the  characteristic  length  should  be  based  on  a 
reference  boundary- layer  thickness.  Hence  we  Introduce  the  following 
scaled  variables: 


<P  -  <t> 

(55a) 

?  =  <p/  <s0 

(55b) 

z  =  z/60 

(55c) 

where 

60  =  »^o/*/,Re 

(56) 

is  a  reference  boundary-layer  thickness 

location  denoted  by  4>0.  Also  we  define 

evaluated  at  the  streamwise 

a  Reynolds  number  based  on  60  , 

R  -  60 Re 

(57) 

In  the  new  scaled  variables,  the  disturbance  equations,  Eqs.  (48)  -  (51), 
without  any  approximation  become 


^Momentum  Equation 

62  l!u  +  ifu  +  h2  -  sn/z  Vo  ~  -  U0R60  ~  +  (mH0  -  l  m$0U0 
3$2  dip2  dz2  3i|>  3<t> 

+  m  *2  V0hK^  -  6jry^)u  -  (m  ^Eo  +  p-  60V0)v 

,  ni  n  2 h  *  i/  h  «  C  P  i~.  ( c  Q  \ 

ijj-Momentum  Equation 

62  iiv  +  ifv  +  h2  Ifv  _  m1/2Vo  av  .  uoR  ^  +  (_mHo  +  1  m6oUo 

3<J)2  dip2  dZ2  dip  dip 

-  m1/2V0h^  -  60  ^)v-  -  (m/z  G0  +  2U0R2hK^)u 

-  m3/z  60  z:  +  2hK,R60  -  h  ^  =  0  (59) 

3l|)  ^  d<P  dip 

z-Momentum  Equation 

6§  “  +  “  +  h2  ^  -  m1/2V0  ~  -  U#R60  —  -  h2  ^  •>  0  (60) 

3(|)2  dp2  dz2  dp  dp  dz 

Continuity  Equation 

R60  —  +  z:  +  h  -  j  60mu  -  hK  v  =  0  (61 ) 

3d)  3d)  3z  w 

where 


u  =  Ui,  v  =  RVj,  w  =  Rwj,  p  =  R2Pi 

(62) 

m  =  <j>o /P 

(63) 

KV  •  5"kg,  “  ■  K*  U 

(64) 

(65) 
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and  Uo»  V0,  E0,  H0,  and  G0  are  functions  related  to  the  basic  flow 
through  Eqs.  (40a)  -  (40e). 

We  identify  60  with  the  pressure  gradient  parameter  defined  by  Eq. 
(33)  and  k^  as  the  curvature  of  the  invlscid  streamlines.  The  value 
of  at  ip  =  0  is  the  curvature  of  the  wall.  In  general,  k^  is  a 
function  of  <j>  and  \jj. 

B.  A  Normal  Mode  Solution 

Equations  (58)  -  (61)  constitute  a  system  of  homogeneous  linear 
partial -differential  equations  with  variable  coefficients.  They  are  of 
the  elliptic  type.  The  class  of  disturbances  considered  here  Is  the 
Gortler  vortices.  That  is,  the  disturbances  are  periodic  in  the  z- 
direction.  The  dependence  on  z  can  be  eliminated  from  the  disturbance 
equations  once  the  wavelength  of  the  periodic  disturbances  is  specified. 
The  result  is  a  system  of  partial-differential  equations  in  and 
as  independent  variables.  Now,  the  boundary  conditions  are 


at 

ip  =  0 

(66) 

as 

\p  -►  cn 

(67) 

Conditions  (66)  are  the  no-slip  and  no- penetration  conditions  for  solid 
walls.  Conditions  (67)  express  the  physical  fact  that  the  disturbances 
must  decay  In  the  free  stream.  For  high  Reynolds  numbers,  the  basic- 
flow  quantities  are  slowly  varying  with  Consequently,  the  basic  flow 
can  be  considered  quasi-parallel,  thereby  allowing  the  separation  of  $ 
from  the  disturbance  equations.  Hence,  we  evaluate  the  basic-flow 
variables  at  a  fixed  location,  say  <j> o,  and  assume  the  solution  in  the 
normal -mode  form 
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(68a) 


u  =  u(^)cosaIe>a  * 

v  =  v^cosaze^^  (68b) 

w  =  w(^)sinaze/ad<i>  (68c) 

p  =  p(^)cosa zefa<®  (68d) 

The  functions  u,  v,  w,  and  p  on  the  right-hand  side  should  not  be 
confused  with  the  total  flow  variables  defined  In  Eqs.  (47a)  -  (47d), 
which  will  not  be  referenced  in  the  rest  of  this  paper.  Substituting 
Eqs.  (68a)  -  (68d)  into  Eqs.  (58)  -  (61),  we  obtain 

j>-Momentum  Equation 

u"  -  V0u 1  +  (H0  -  h2oc2  -  alio  -  \  Mo)u  -  E„v  +  rx  ■  0  (69) 

ip-Momentum  Equation 

v"  -  V0V  +  (-H0  -  hV  -  aU0  +  \  6oU0)v  -  (2U0hG^  +  G0)u 

-  BoU1  -  hp1  +  r2  =  0  (70) 

z-Momentum  Equation 

w "  -  V0w'  +  (-n2a2  -  oU0)w  +  h2ap  +  r3  =  0  (71) 

Continuity  Equation 


v'  +  (a  -  j  Bo)u  +  haw  +  l\  =>  0 


where 


a  =  Ra 


r  _  /u  ui/  t  in  do  A  a2  6oV2h\ 
Fl  -  (V°hVR  ~  +  R ^--JLh~)u 


-  It  (BoVo  +  2hK  a)v  +  If  v' 
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(75) 


«  ( -  V  o  h  K  +  4^  —  +  £ 
V  K  dp  K 


6o  -r^-)v  +  2hK  au 


<4*  - +  £  )« 

K  d<t>  K 


r4  =  -  hK, v 


Gw  =  R2K, 

N  \p 


and  prime  denotes  differentiation  witn  respect  to  j7. 

The  variation  of  the  streamline  curvature  in  the  ^-direction  (i.e., 
normal  to  the  wall)  is  included  in  the  function  K^.  Hence,  in  general, 

Gjjj  is  a  function  of  tp.  We  identify  the  value  of  at  ip  -  0  as  the 

Gortler  Number.  Recalling  the  definitions  of  R  and  K^,  Eqs.  (57)  and 
(64),  we  see  that  the  present  definition  of  Gortler  number  is  based  on 
the  reference  boundary- layer  thickness  60,  while  the  original  parameter 
identified  by  Gortler  is  based  on  the  momentum  thickness.  However,  the 
two  thicknesses  are  related  (e.g.,  for  the  Blasius  profile  6m  *  0.664  60). 

The  boundary  conditions,  Eqs.  (66)  and  (67),  become 

u=v=w=0  at  p  =  0  (79) 

u,  v,  and  w  -*■  0  as  ij>  -*•  »  (80) 

The  homogeneous  systen  defined  by  Eqs.  (69)  -  (72),  (79),  and 
(80)  is  an  eigenvalue  problem.  The  characteristic  equation  of  this 
eigenvalue  problem  may  have  the  form 

T(a,  o,  R)  =  0  (81) 


where  a  is  the  wavenumber,  a  i?  the  growth  (amplification)  rate,  and  R 
is  a  reference  Reynolds  numbet  .  Because  of  the  dependence  of  the  basic 
flow  on  the  wall  geometry  and  the  explicit  existence  of  the  metric 
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coefficient  h  and  its  derivatives  in  Eqs.  (69)  -  (72),  the  surface  given 
by  Eq.  (81)  Is  configuration  dependent.  That  is.  In  general,  each 
geometry  may  have  different  characteristic  surfaces  In  the  space  o,  a,  and 
R. 

C.  An  Equivalent  System  of  First-Order  Equations 

Equations  (69)  -  (72)  are  reduced  to  a  system  of  first-order 
differential  equations  by  letting 

Zi  =  u,  z2  =  u',  z3  =  v,  z4  =  w,  z5  »  w'  and  z6  =  hp 


The  result  is 


z'  =  Az 

where  z  is  the  vector 

2  s  (Zi ,  z2,  z3,  Z4,  Z5,  Zs t ) 


(82) 


(83) 


and  A  is  a  matrix  that  can  be  written  as  the  sum  of  two  matrices  according 
to 


A  =  A0  +  AA 


a2+aUo-Ho+  \  6oU0 
1 


An  = 


2  D0 


-a+  o  6 
0 
0 
-A 


1 

Vo 

0 

0 

0 

1 


0 

Eo 

0 

0 

0 


0 

0 

-a 

0 

a2+aUo 


-0-  2  80  -(a2+aUo+H0-  j  8oU0)  0V0 


(84) 

0  0 

0  0 

0  0 

1  0 

Vo  -a 

-a  0 


where 


A  =  2U0G2 


”  aV0  Sq  ^  o  60V0 


(85) 
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and  the  matrix  AA  is  given  in  Appendix  C.  We  note  that  the  matrix  A0 
contains  all  the  terms  of  0(1),  while  the  matrix  aA  contains  all  the 
higher-order  terms.  The  boundary  conditions  for  the  system  (82)  are 


Zi  =  z3  =  z4  =  0 

at 

ip  -  0 

(86) 

Zi ,  z3 ,  and  z4  •*  0 

as 

Ip  -*•  CO 

(87) 

We  recall  the  fact  that  the  Gortler  number  in  Eq.  (85)  is  a  function 
of  i j>. 

The  above  reduction  of  the  governing  equations  to  a  system  of  six 

first-order  equations  requires  the  first  derivatives  of  the  basic-flow 

12  13 

quantities.  In  the  work  of  Smith  and  Floryan  and  Saric  ,  the  dis¬ 
turbance  equations  were  reduced  to  two  coupled  equations,  one  of  second 
order  in  the  u-component  and  the  other  of  fourth  order  in  the  v-compo- 
nent.  Then,  Floryan  and  Saric  reduced  the  new  equations  to  a  system  of 
six  first-order  equations.  However,  their  method  requires  some  of  the 
second-order  derivatives  of  the  basic  Flow,  which  demands  a  higher 
accuracy  of  the  basic-flow  solution  than  it  is  required  in  the  present 
method. 

We  define  a  basic  system  as 

z'  -  A0z  (88) 

subject  to  the  boundary  conditions  (86)  and  (871.  We  note  that  the 

13 

zeroth-order  disturbance  equations  developed  by  Floryan  and  Saric 
for  the  case  of  a  circular  arc  can  be  put  in  the  form  (88).  A  method 
for  solving  the  basic  system  is  presented  in  Section  IV.  However,  the 
eigenvalues  of  the  basic  system  will  be  slightly  different  from  those  of 
the  original  system,  Eqs.  (82),  because  the  matrix  AA  was  not  Included. 
In  Section  V,  a  perturbation  procedure  is  presented  to  obtain  the 
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correction  to  the  eigenvalues  due  to  a  small  deviation  of  the  coefficient 
matrix  from  a  basic  value. 

IV.  A  Method  For  Solving  The  Basic  System 

The  system  (88)  subject  to  boundary  conditions  (86)  and  (87)  is  a 
two-point  boundary-value  problem.  As  mentioned  previously,  this  is  an 
eigenvalue  problem.  Any  one  of  the  parameters  a,  a,  and  R  can  be 
treated  as  the  eigenvalue  when  the  other  two  parameters  are  prescribed. 

An  alternative  to  the  parameter  R  is  the  Gortler  number  G^(0)  based  on 
the  wall  curvature,  in  which  case  we  rewrite  Eq.  (85)  as 

A  =  2U0Gjj(0)rffi)  -  oV0  +  Go  +  \  Mo  (89) 

where  r(0)  =1.  If  the  deviation  of  the  curvature  of  the  streamlines 
from  the  value  at  the  wall  is  neglected,  then  r  =  1  for  all  ip  and  G^(0) 
is  the  eigenvalue.  However,  for  the  purpose  of  the  present  section,  r 
is  assumed  to  be 

r  =  r(lp)  for  0  <  ip  < 

(90) 

r  =  rffim)  ip  '> 

The  difference  between  f(ij;)  and  r(om)  tor  .  ,  will  oe  accounted  for 

by  the  perturbation  technique  presented  in  Section  V. 

A .  Boundary  Conditions  at  I nfjnity 

Application  of  the  asymptotic  boundary  condition  (87)  at  a  finite 
value  of  Ip  may  lead  tc  inaccurate  results.  However,  outside  the  boundary 
layer  (i.e.,  \p  >  10),  the  system  of  eouations  has  constant  coefficients 
because 
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Uo  1 ,  H0  ■+  0,  V o  Vo ,  E0  0, 

and  Go  •*  Go  as  ip  -*  °°, 

where  V0  and  G0  are  constants.  Thus,  we  write  the  system  of  constant 
coefficients  as 

2 1  =  Cz  (91) 


where 


C  = 


0 

,  1 

a2  +  a  +  j  Bo 


-a 

0 

0 

-A 


and 


0 

Vo 

0 

0 

0 


0 

0 

0 

0 

0 


-a  -  o-  Bo  -(a2  +  a  -  p-  Bo) 


0 

C 

-a 

0 

a2  +  a 
a 


1  1/ 


0 

0 

0 

1 

Vo 

-a 


0 

0 

0 

0 

-a 


°J 


A  =  2G2(0)r(^m)  -  aV0  +  G0  +  j  BoV0 


(92) 


(93) 


Let  us  denote  the  characteristic  roots  of  C  by  Ai,  A2,  ....  A6 
where  the  An  are  given  by 

A2  -  V0A  -  (a2  +  a  +  Bo)  =  0  (94a) 

A1*  -  VoA3  -  (2a2  +  a)A2  +  a2VoA  +  a2(a2  +  a  -  ^  So)  3  0  (94b) 

We  consider  the  case  in  which  three  of  the  Ap  are  positive  real  numbers, 
while  the  other  three  are  negative  real  numbers.  Let  Ai,  A2,  and  A3 
be  the  positive  roots.  When  Bo  =  0,  the  characteristic  roots  of  C 
reduce  to 

a,  A2,  A3,  -  a,  A5,  and  A6  (95) 
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where 


X2  =  X3  3  J  [Vo  +  /V§  +  4(a 2  +  or)]  (96a) 

and 

X5  =  X6  =  -  \  [-Vo  +  /v?  +  4(o2  +o)]  (96b) 

For  the  special  case  a  =  aV0,  Eq.  (96b)  gives 

X5  =  X6  =  “  d 

so  that  -  a  is  a  root  of  multiplicity  3.  When  a  <  -  ^-(Vo  +  4a2) 

Eqs.  (96)  show  that  X2,  X3,  X5,  and  X6  are  complex.  This  case  of 
complex  roots  is  not  treated  in  the  present  work  as  mentioned  earlier. 

Now,  an  analytic  solution  to  the  system  (91)  Is  easily  obtained  in 
terms  of  elementary  functions.  The  system  (89)  has  six  linearly  indepen¬ 
dent  solutions.  The  general  solution  Is  obtained  by  superposition  of 
these  fundamental  solutions.  Hence,  we  write 

z  =  Yb  (97) 

where  the  ith  component  of  the  vector  b  is 

bi  =  ci  exp(X14)  (98) 

The  fundamental  matrix  Y  is  not  a  constant  matrix  in  the  case  of  repeated 
roots.  It  includes  linear  terms  in  ^  when  60  =  0;  Tor  the  special  case 
when  -  a  has  a  multiplicity  of  3,  it  includes  quadratic  terms  In 
Therefore,  its  value  depends  on  ijJ. 

Because  of  the  linear  independence  of  the  fundamental  solutions, 
the  system  (97)  can  be  inverted  to  yield 

Y-1z  =  b  (99) 

Now  the  asymptotic  boundary  condition  (87)  requires  that  bi  ■  b2  3  b3  *  0, 
as  can  be  seen  from  Eq.  (97)  and  (98).  It  follows  from  Eq.  (99)  that 


/ 
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Tz  =  {0}  (100) 

where  T  is  a  3  x  6  rectangular  matrix  given  by 

Ti,r[v'l]u  j :  1: 2:.?6  <,o,> 

where  [Y-1]..  denotes  the  element  ij  of  the  matrix  Y"1. 

*  J 

By  evaluating  T  at  a  finite  value  of  ij)  -  ip  ,  Eq.  (100)  provides 
three  boundary  conditions  in  the  form  of  linear  combinations  of  the 
values  of  the  variables  zi,  . . . ,  z6  at  i|>  =  ij;  .  This  form  suits  the 
method  used  to  integrate  the  basic  system  for  ip  < 

We  note  that  the  matrix  T  in  Eqs.  (100)  depends  on  ip_  when  Bo  =  0, 
which  is  calculated  by  inverting  a  6  x  6  matrix  for  each  iteration  and 
for  each  eigenvalue.  Next,  we  develop  an  alternative  form  for  T  that  is 
independent  of  ip.  To  accomplish  this,  we  reduce  the  matrix  C  to  a 


Jordan  canonical  form  by  introducing  the  similarity  transformation 


J  =  P"*CP 


(10?) 


When  Bo  t  0,  the  Xn  are  distinct  and  their  corresponding  eigenvectors 
are  orthogonal.  In  this  case,  J  is  a  diagonal  matrix  and  P  is  a  matrix 
whose  columns  are  the  eigenvectors  of  C.  We  arrange  these  columns  so 
that  the  positive  eigenvalues  appear  in  the  first  three  rows.  When 
Bo  =  0,  X2  =  X3  and  X5  =  X6,  and  hence  J  has  the  form 


X2  T 1  0 


0  0 


0  0 
0  0 


0  0 


X5  T2 


(103) 
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where  ti  and  t2  are  constants  that  we  do  not  need  to  specify.  In  this 
case,  P  is  defined  in  Appendix  D.  Letting 


z  =  PE,  (104) 

in  (91)  and  using  (102),  we  obtain 

V  =  JC  (105) 

We  note  that  the  characteristic  roots  in  J  are  arranged  such  that  the 
first  three  roots  are  the  positive  roots.  Therefore,  for  decaying 
solutions  as  ip  ■*  °°,  we  must  impose  the  conditions 

Ci  =  C2  =  Cs  =  0  (106) 

Since  5  =  P~‘z  according  to  Eq.  (104),  conditions  (106)  yield 


Tz  =  {0} 


(107) 


where  T  consists  of  the  first  three  rows  of  P”1. 

Although  the  matrix  T  of  Eq.  (107)  does  not  depend  on  ip.  it  is 
calculated  by  inverting  the  6x6  matrix  P  for  each  iteration  and  for 
each  eigenvalue.  Next,  we  use  the  concept  of  adjoint  to  develop  a 
matrix  T  that  does  not  require  inversion  of  a  matrix.  To  accomplish 

*t  _ 

this,  we  multiply  Eq.  (88)  with  z  (\p) ,  Integrate  the  result  by  parts 
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from  ijj  =  0  to  and  obtain 


(^—  +  z*TAo  )zd\p  =  0 

dip 


(108) 


The  equations  defining  the  adjoint  z*  are  obtained  by  setting  the  coef¬ 
ficient  of  z  in  the  integrand  in  Eq.  (108)  equal  to  zero.  The  result  is 


which  upon  transposition  becomes 


dz 

dljj 

Then,  Eq. 


T  * 

+  Ajz  -  0 
(108)  becomes 


(109) 


(110) 


K  io 

We  recall  that  the  boundary  conditions  at  -  0  are  defined  in  Eq. 
(86),  while  the  boundary  conditions  at  are  defined  in  Eq.  (107), 
where  T  consists  of  the  first  three  rows  of  P  ’ .  We  choose  the  adjoint 
boundary  conditions  such  that  the  two  terms  in  Eq.  (110)  vanish 
independently;  that  is 


*Tz  -  0 

at 

*  =  *m 

(111) 

*Tz  =  0 

at 

ip  =  0 

(112) 

In  expanded  form,  Eq.  (112)  becomes 


^  if  if  ic  *  *  _ 

ZiZi  +  z2z2+  Z3Z3  +  Z4Z4  +  Z5Z5  +  z6z6  =  0  at  \p  =  0 


which  upon  using  Eq.  (86)  reduces  to 

if  if  ic  _ 

z2z2  +  Z5Z5  +  z6z6  =  0  at  4>  = 


0 


(113) 


We  choose  the  boundary  conditions  at  ip  =  0  such  that  each  of  the  coef¬ 
ficients  of  z2(0),  z5(0),  and  z6(0)  vanish  independently.  The  result  is 

z2  =  z*  =  z*  =  0  at  ip  =  0  (114) 

To  define  the  boundary  conditions  at  ip  =  i pm,  we  view  Eq.  (104)  as  a 
linear  transformation  relating  C  and  z  at  ip  =  ipm.  Substituting  Eq. 

(104)  Into  (111),  we  have 
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z*TP£  =  0  at  ?  =  4>m 

or 

(PTz*)T5  =  0  at  *  =  (115) 

Letting 

★  T  * 

t  =  P  z  (116) 

we  rewrite  Eq.  (115)  as 

C*TS  -  0  at  y  =  (117) 

or 

£l£l  +  &2C2  +  ^3^3  +  £■.£■.  +  S5C5  +  ^6^6  =  0  at 

(118) 

Using  the  boundary  conditions  (106),  which  are  equivalent  to  the  boundary 
conditions  (107),  we  reduce  Eq.  (118)  to 

+  CsZs  +  tUe  =0  at  ?  «  (119) 

Since  £4,  ?5,  and  are  linearly  independent  because  P  is  a  nonsingular 
matrix,  each  of  their  coefficients  must  vanish  independently;  that  is 

★  ★  ★  _  _ 

=  0  at  <|i  =  ^  (120) 

Then,  it  follows  from  Eq.  (116)  and  (120)  that  the  adjoint  boundary 
conditions  at  ip  =  IjJ  are 

T*z*  =  0  at  ^  (121  ) 

★  T  T 

where  T  consists  of  the  last  three  rows  of  P  when  60  =  0  and  Y  when 

80  t  0.  We  note  that  the  calculation  of  T  does  not  involve  any  matrix 

inversion.  Thus,  the  adjoint  problem  is  defined  by  Eq.  (109)  subject 

to  the  boundary  conditions  (114)  and  (121). 
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To  determine  boundary  conditions  for  z  that  do  not  Involve  inver¬ 
sion  of  matrices,  we  consider  the  adjoint  of  the  adjoint  problem.  When 
ip  >  "i£m,  Eq.  (109)  becomes 

+  C‘z  =0  (122) 

d^) 

*  T 

As  before,  we  reduce  the  matrix  C  =  -  C  into  a  Jordan  canonical  form 
by  using  the  similarity  transformation 

J  =  P  ’CP  (123) 


When  3o  t  0,  the  An  are  distinct  and  their  corresponding  eigenvectors 
are  orthogonal.  In  this  case,  J*  is  a  diagonal  matrix  and  P*  is  a 
matrix  whose  columns  are  the  eigenvectors  of  C*.  We  arrange  these 
columns  so  that  the  positive  eigenvalues  appear  in  the  first  three  rows. 
When  0o  =  0,  J*  has  the  form 


(124) 


★  ★  ★  * 
where  ii  and  x2  are  constants  and  A2  =  -A5  and  A5  =  -A2.  In  this  case, 

P*  is  defined  in  Appendix  D.  Then,  using  the  transformation 
★  ★  * 

z  =  P  5  (125) 


we  rewrite  (122)  as 

^  =  jY 

dip 

Hence,  for  decaying  solutions,  we  must  put 


(126) 
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5?  =  S*  =  et  =  0  at  ip  =  <J»m  (127) 

Then,  following  steps  similar  to  those  in  Eq.  (108)  -  (121),  we  find 
that  the  boundary  conditions  for  z  at  ip  =  <|>m  are 

Tz  =  0  at  ^  (128) 

*T 

where  T  consists  of  the  last  three  rows  of  P  ,  which  do  not  involve 
inversion  of  matrices. 

B.  Boundary  Conditions  at  tne  Mall 

The  boundary  conditions  (86)  and  (114)  can  be  rewritten  in  matrix 
form  as 

Wz  =  {0}  at  ^  =  0  (129) 

for  the  basic  system,  and 

W*z  =  {0}  at  =  0  (130) 

★ 

for  the  adjoint  system.  Here,  W  and  W  are  3x6  rectangular  matrices 
with  all  elements  being  zero  except  that 

W i i  =  W23  =  W34  =  W12  =  W25  -  W36  =  1 

In  summary,  the  basic  system  is  defined  by  Eqs.  (88),  (128),  and  (129), 
and  the  adjoint  system  is  defined  by  Eqs.  (109),  (121),  and  (130). 

C.  Numerical  Procedure 

Numerical  integration  of  these  two-point  boundary-value  problems 
are  obtained  by  using  a  computer  code  (SUP0RT)  developed  by  Scott  and 
Watts  .  The  method  of  integration  uses  superposition  coupled  with  an 
orthononmalization  procedure  and  a  variable-step  Runge-Kutta-Fehlberg 
scheme.  We  fix  the  values  of  any  two  of  the  parameter  a,  a,  and  G^,  and 
consider  the  third  one  as  an  eigenvalue.  We  make  an  initial  guess  for 


33 


the  eigenvalue.  Starting  at  the  position  we  carry  out  the 

integration  toward  the  wall  <p  =  0.  One  of  the  boundary  conditions  at 
the  wall  will  not  be  satisfied  if  the  guessed  eigenvalue  is  not  the 
exact  value.  A  Newton-Raphson  iteration  process  is  used  to  obtain  a 
better  approximation  of  the  eigenvalue.  The  computer  code  is  designed 
such  that  the  boundary  condition  that  has  to  be  Iterated  on  is  the  last 
specified  one.  Hence,  for  the  basic  system  we  iterate  on  the  condition 
z4  =  0  at  the  wall.  The  condition  of  convergence  was  taken  as 

|z4(0)/z4  |  <  10"6  (131) 

A 

where  z4  is  the  maximum  value  of  z4  in  the  domain  of  numerical  integra¬ 
tion.  For  the  adjoint  problem  the  condition  was  taken  as 

|zl«»/i:  |  <  10'6  (132) 

Conditions  (131)  and  (132)  were  used  for  all  the  calculations  presented 
in  this  paper.  However,  for  some  of  the  cases,  higher  accuracy  was  used 
in  Eqs .  (131)  and  (132).  However,  the  change  in  the  eigenvalues  and 
eigenfunctions  are  of  0( 1 0~ 6 ) . 

V.  A  Perturbation  Procedure  to  Account  for  Small  Deviations  in  the 

Matrix  A0 

The  matrix  A  in  Eq.  (84)  is  expressed  as  the  sum  of  two  matrices 
A0  and  AA,  where  the  elements  of  AA  are  small  compared  with  those  of  A0. 
In  Section  IV  we  presented  a  method  of  solution  for  the  basic  system 
whose  matrix  Is  A0.  New,  we  present  a  perturbation  approach  to  find  the 
change  in  the  eigenvalue  due  to  a  small  change  Ai  in  A0.  However,  we 
need  not  identify  the  change  Ai  as  the  matrix  AA  in  Eq.  (84).  The 
matrix  A,  represents  a  general,  however,  small  change  In  A0;  Aj  may 
incorporate,  in  addition  to  AA,  changes  in  the  basic-flow  quantities  due 
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to  displacement  speed  and/or  viscous  curvature  effects.  Since  the 
perturbation  in  the  eigenvalue  due  to  each  higher-order  effect  is 
small,  then  interaction  between  these  effects  in  negligible.  Hence,  one 
may  study  the  influence  of  each  higher-order  term  independent  of  the 
others.  This  will  point  out  which  second-order  terms  are  the  most 
stabilizing  or  destabilizing.  This  is  helpful  in  the  aerodynamic 
design  of  surfaces  with  regions  of  concave  curvature. 

We  consider  the  eigenvalue  problem 

z‘  =  (A0  +  eAi )z  (133) 

with  the  boundary  conditions  (86)  and  (87).  We  assume  that  e  «  1. 

First  we  solve  the  basic  system 

zj  =  A0z0  (134) 

with  the  same  boundary  conditions.  Let  X0  denote  the  eigenvalue  of  this 
system.  We  expand  the  solution  for  system  (133)  in  the  form 

z  =  z0  +  ec  (135) 

and  the  eigenvalue  as 

X  =  Ao  +  eXj  (136) 


Since  A0  depends  on  the  eigenvalue  X,  we  use  the  Taylor-series  expansion 

9  A  o 


A0 (Ao  +  eXi )  -  Ao (X0 )  + 


ax 


eXi  +  ... 


037) 


Substituting  expansions  (135)  and  (137)  into  Eq.  (133),  neglecting  terms 
0(e2),  and  using  Eq.  (134),  we  obtain 

3A0 


C1  -  A0c  =  (Ai 


ax 


+  A,)z0 


(138) 


The  boundary  conditions  (86)  and  (87)  give 


Cl  =  C3  =  C-  =  0 

at 

if  =  0 

(139) 

Cl  *  C3  t  C4  •*  0 

as 

\p  -*•  00 

(140) 

Since  the  homogeneous  problem  obtained  from  the  system  (138)  -  (140) 

has  a  nontrivial  solution,  namely  z0»  then  the  nonhomogeneous  problem 
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has  a  solution  if  and  only  if  a  solvability  condition  is  satisfied  . 
To  obtain  this  condition,  we  multiply  Eq.  (138)  by  the  solution  of 
the  adjoint  homogeneous  problem  and  integrate  the  result  ^rom  if  -  0  to 
if  ■*  The  result  is 


(141) 


Hence,  the  eigenvalue  of  the  original  system  (133)  becomes 


(142) 


The  integrals  in  Eq.  (142)  have  to  be  evaluated  numerically.  In  the 
interval  0  <  if  <_  ifm  the  values  of  z0  and  z0*  are  known  from  the 
numerical  solution  of  the  basic  problem  and  its  adjoint.  For  if  >  If 
the  analytic  solution  outside  the  boundary  layer  is  used.  The  eigen¬ 
value  X  in  Eq.  .142)  can  be  identified  as  a,  a,  or  R  and  its  alterna¬ 
tive  G^(0).  However,  from  the  structure  of  the  matrix  A0»  it  is  advan¬ 
tageous  to  identify  X  with  G^(0)  since  in  this  case  all  elements  in  the 
matrix  9A0/9X  are  zero  except  one  element.  Let  denote  an  element 


1 
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(143) 


of  3Ao/9A,  hence  D,,  =  0  for  all  i  and  j  except 

J 

D6 i  =  -  2U0r 

where  r  is  defined  by  Eq.  (90).  Thus,  Eq.  (142)  gives 


An  + 


a 

{ 


Zq  eAlZodil/ 


1. 


(144) 


rUoZ6Zidip 


where  Zi  is  the  first  component  of  z0  and  z6  is  the  sixth  component  of 
Zo-  In  Eq.  (144),  A  and  A0  stand  for  Gj^(O)  with  and  without  the 
correction  due  eAi,  respectively. 

As  a  first  application,  we  consider  the  correction  due  to  restricting 
the  function  r(^)  to  be  a  constant  for  ijj  >  Let  B^.  denote  an  element 

of  eAi,  then  B..  =  0  for  all  i  and  j  except 

'  J 


>6  1 


{ 


for  ip  < 

—  m 


(145) 


2U0A0 [r(i^m)  -  r(ij»)]  for  \p  > 


We  note  that  l)0  =  1  for  \p  >  \jj  .  Equation  (144)  gives 


A/Ao  =  (Ii  +  1 3 )  /  ( 1 1  +  Iz) 


where 


r  . 

Ii  =  J  rUoZeZitty 

0 


★  _ 

rz6Zid^ 


(146) 


(147a) 


(147b) 


00 

J 


Is  =  r(i|»m)  J  ZsZidip 


(147c) 
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As  a  second  application  of  Eq.  (144),  we  consdier  the  effect  of 
displacement  thickness  of  the  Blasius  flow.  It  follows  from  Eqs.  (45)  and 
(46)  that  the  perturbations  in  V0  and  G0  are 


Vi 

=  -  V0[l  -  Real ( t ‘ 7 2 ) ] 

(148) 

G, 

=  1  V0[l  -  Real ( t 3 7 2 ) 3 

(149) 

where 

t  = 

:  (1  +  iif)//Re<J>0  )  1 

(150) 

The  nonzero  elements  of  eAi  are  B22  =  V2,  B55  =  Vlf  B61  =  aV i 

-  Gi , 

and  B64 

=  aVi.  Then  Eq.  (144)  gives 

A 

/,C°  ★  ★  ★  * 

=  A0  +  ]  [  V  i  (z2z2  +  z5z5  +  az4z6  +  aZjzJ 

w  0 

CO 

-  G i z i z & ] J*  2UoZiZ6di); 

(151) 

The  destabilizing  effect  of  Gi  is  clear  from  Eq.  (151)  since  Its 
contribution  to  the  integral  is  negative. 

VI.  Results  and  Discussion 

Gortler^  showed  that  the  centrifugal  instability  of  boundary-layer 
flows  along  concavely  curved  walls  is  primarily  a  function  of  the  momen¬ 
tum  thickness  of  the  basic-flow  profile  and  not  of  the  details  of  the 
shape  of  the  boundary  layer.  We  recall  that  in  Gortler's  model  the 

transverse  component  of  the  basic  flow  was  neglected.  This  component 

12  13 

was  accounted  for  in  the  models  of  Smith  and  Floryan  and  Saric  .  Very 

accurate  numerical  calculations  of  the  neutral  stability  curves  of  these 

1 3 

two  models  for  the  Blasius  profile  were  provided  by  Floryan  and  Saric  . 
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Their  results  demonstrated  the  destabilizing  effect  of  the  transverse 
component  of  the  basic  flow,  especially  for  small  wavenumbers. 

A.  The  Blasius  Flow 

Figure  4  shows  the  neutral  stability  curves  for  the  Blasius  flow 

obtained  for  different  models.  Curve  1  in  Fig.  4  is  Gortler's  model. 

The  pertinent  equations  are  obtained  from  the  basic  system,  Eq.  (88), 

by  setting  Vo  =  Go  =  H0  =  0  for  all  ip  ( i . e. ,  we  neglect  the  terms  due  to 

the  transverse  component  of  the  basic  flow).  Curve  2  is  the  present 

1 3 

model,  Curve  3  is  the  Floryan-Saric  model,  and  Curve  4  is  a  modified 
Smith  model.  The  governing  equations  of  these  three  models  are  identi¬ 
cal,  Eq.  (88).  The  differences  among  the  neutral  stability  curves  of 
these  models,  as  shown  in  Fig.  4,  are  due  to  different  treatments  of  the 
term  Go-  In  the  present  model.  Curve  2,  G0  -  -1/2  V0  3  -0.4302  as 
ip  -*  as  can  be  seen  from  Eq.  (41  e).  Outside  the  boundary  layer 

(i.e.,  u  C  1.0),  the  terms  V0  and  G0  represent  the  slope  and  curvature 
of  the  viscous  streamlines,  respectively.  As  discussed  in  Section  II, 
these  terms  do  not  vanish  outside  the  boundary  layer  in  both  the  first- 
and  second-order  boundary-layer  solutions.  However,  as  ri  »  they 

vanish  only  in  the  second-order  solution  when  we  account  for  the  dis- 

13 

placement-thickness  effect.  Floryan  and  Saric  ,  in  the  calculation  of 

Curve  3,  put  3V0/34>i  3  0  for  i p  >  ipm,  where  cpi  =  cv<p  and  cv  =  1/R. 

Equation  (39e)  shows  that  the  part  that  does  not  vanish  as  n  -►  00  in  the 

expression  of  3hv/34>  is  the  first  term,  namely  (V  +  nF).  When  we 

neglected  this  term  everywhere  (i.e.,  for  all  i|/),  we  obtained  a  neutral 

stability  curve  that  is  indistinguishable  from  Curve  3  of  Floryan  and 

1 2 

Saric.  We  recall  that  in  the  Smith  model,  the  term  9hv/34>  was  not 
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included  in  the  analysis.  Hence,  in  the  modified  Smith  model,  Curve  4, 
we  let  G0  =  0  for  all  4). 

To  explain  the  effect  of  Ga  on  the  neutral  stability  curves,  we 
recall  Eqs.  (89)  and  (96),  rewritten  here  for  r  =  1  and  a  =  0, 

A  =  2U0G^  +  Go  (152) 

X2  =  -  \  C-Vo  +  /Vl  +  4a"2]  (153) 

where  V0  =  0.8604  for  the  Blasius  flow.  The  disturbances  aecay  like 
exp(-|X2|^)  outside  the  boundary  layer.  From  Eq.  (153),  we  see  that  for 
small  wavenumbers,  a  <  0.1,  |X2|  «  1,  and  hence  the  disturbances  decay 
very  slowly  as  4>  °°.  Therefore,  changes  in  the  basic-flow  quantities 

outside  the  boundary  layer  affect  the  neutral  stability  characteristics 
for  small  wavenumbers.  Furthermore,  the  term  G^  is  proportional  to  the 
curvature  of  the  inviscid  streamlines.  The  presence  of  G0  with  G^  In 
the  same  term,  attaches  the  meaning  of  curvature  to  G0.  Therefore,  we 
can  view  A  in  Eq.  (152)  as  a  measure  of  the  effective  curvature  that 
includes  viscous  effects.  We  note  that  for  the  Blasius  flow  and  non¬ 
separating  Falkner-Skan  flows,  C0  is  negative  for  all  i p,  as  can  be  seen 
from  Eq.  (40e).  It  follows  from  Eq.  (152)  that  G^  takes  a  larger  value 
when  we  include  the  negative  quantity  G0,  Curve  2,  than  when  we  neglect 
it  completely,  as  in  the  modified  Smith  model,  Curve  4.  For  Curve  3,  G0 
is  slightly  negative  and  it  vanishes  outside  the  boundary  layer.  How¬ 
ever,  for  large  wavenumbers  or  growing  vortices  (i.e.,  a  >  0)  Eqs, 

(153)  and  (96)  indicate  that  the  disturbances  decay  very  fast  outside 
the  boundary  layer.  Hence  the  effect  of  the  term  G0»  which  has  more 
influence  outside  the  boundary  layer  than  inside  it,  should  diminish. 
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This  also  follows  from  Eqs.  (152)  when  the  Gortler  number  GN  takes  large 
values. 

Finally,  we  note  that  =  12  was  used  for  all  the  calculations 
presented  in  Fig.  4.  The  accuracy  of  the  solution  was  tested  against 
changes  in  ij>  .  Increasing  ^  up  to  24  produced  insignificant  changes, 
0(1  Cf6),  in  both  the  eigenvalues  and  the  eigenfunctions.  Also,  the 
effect  of  the  accuracy  of  the  basic-flow  solution  was  evaluated  by  using 
a  second-order  finite-difference  scheme  with  An  =  0.01  Instead  of  the 
fourth-order  Runge-Kutta  scheme  with  An  =  0.01.  For  the  test  case 
a  =  0.01,  the  change  in  the  eigenvalue  was  insignificant  0(10-5),  while 
the  changes  in  the  maximum  values  of  the  eigenfunctions  were  about  7%. 
This  is  consistent  with  the  fact  that  the  eigenvalue  is  a  global  value 
while  the  eigenfunction  is  a  point  solution.  As  a  further  test  on  the 
eigenvalues  obtained  for  all  the  cases  presented  in  this  chapter,  the 
adjoint  problem  was  solved  and  the  obtained  eigenvalues  were  compared 
with  those  of  the  basic  system.  The  discrepancy  between  the  two 
eigenvalues  thus  obtained  was  0(10  5). 

Table  1  gives  the  Gbrtler  number,  G^,  as  a  function  of  the  wave- 
number,  a,  for  the  neutral  stability,  a  =  0,  of  the  Blaslus  flow  for  the 
previously  mentioned  models. 

B.  Effect  of  Displacement  Thickness 

The  perturbation  procedure  presented  In  Section  V  was  used  to 
account  for  the  effect  of  changes  in  the  basic  flow  due  to  the  displaced 
body  on  the  neutral  stability  of  the  Blasius  flow.  Eouations  (148)- 
(150)  show  that  the  corrections  Vi  and  Gi  depend  explicitly  on  the 
Reynolds  number  Re.  We  let  Rex  denote  Re<t>0.  Figure  5  shows  the  neutral 
stability  curves  (the  dashed  lines)  for  the  Blasius  flow  with  correction 
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Table  1 


Values 

Neutral 

of  Gortler 

Stability  of 
Number  G^  as 

the  Blasius  Flow, 

a  Function  of  Wavenumber  a 

gn 

a 

Gortler 

Present 

Floryan- 

Modified 

Model 

Model 

Saric 

Model 

Smith 

Model 

1.00 

3.0392 

2.8238 

2.7954 

2.7629 

0.90 

2.6492 

2.4366 

2.4028 

2.3662 

0.80 

2.2919 

2.0829 

2.0423 

2.0008 

0.70 

1.9669 

1.7623 

1.7175 

i.6t:r 

0.60 

1.6735 

1.4747 

1.4136 

1 . 3605 

0.50 

1.4110 

1.2200 

1.1429 

1.0834 

0.40 

1.1788 

0.9981 

0.8990 

0.8332 

0.30 

0.9761 

0.8100 

0.6792 

0.6088 

0.25 

0.8856 

0.7290 

0.5773 

0.5060 

0.20 

0.8024 

0.5672 

0.4797 

0.4090 

0.15 

0.7264 

0.5947 

0.3853 

0.3176 

0.10 

0.6577 

0.5418 

0.2912 

0.2305 

0.09 

0.6449 

0.5323 

0.2720 

0.2133 

0.08 

0.6324 

0.5233 

0.2525 

0.1961 

0.07 

0.6202 

0.5146 

0.2325 

0.1788 

0.06 

0.6083 

0.5063 

0.2119 

0.1613 

0.05 

0.5967 

0.4983 

0.1905 

0.1435 

0.04 

0.5855 

0.4907 

0.1677 

0.1249 

0.03 

0.5746 

0.4835 

0.1430 

0.1053 

0.025 

0.5693 

0.4800 

0.1295 

0.09482 

0.020 

0.5641 

0.4766 

0.1149 

0.08365 

0.015 

0.5590 

0.4732 

0.09877 

0.07144 

0.010 

0.5539 

0.4700 

0.08001 

0.5752 
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due  to  displacement  speed  for*  Rex  =  103  and  104.  It  Is  clear  from  Fig. 

5  that,  for  a  given  Reynolds  number  and  wavenumber,  the  displacement- 
thickness  effect  is  destabilizing  for  the  Blaslus  flow.  This  effect  Is 
pronounced  for  small  wavenumbers,  which  is  expected  because  the  distur¬ 
bances  extend  farther  outside  the  boundary  layer,  and  there  the  changes 
in  the  basic  flow  due  to  the  displacement  body  are  important.  Figure  6 
shows  the  ratio  of  the  Gortler  number  with  displacement  thickness  effect 
accounted  for  to  the  Gortler  number  when  this  effect  is  absent.  For  a 
given  Reynolds  number,  this  ratio  can  be  viewed  as  the  ratio  of  the 
square  root  of  the  corresponding  critical  curvatures. 

C.  Effect  of  Decaying  Streamline  Curvature 

Herbert11  presented  results  for  the  Gortler  model  with  streamline 
curvatures  decaying  exponentially  away  from  the  wall.  Therefore,  the 
function  r(«f7)  in  Eq.  (89)  has  the  form 


where  c  is  a  constant.  In  this  case,  an  analytic  solution  for  the 
system  (91)  can  be  obtained.  However,  in  a  more  general  situation,  the 
variation  of  the  streamline  curvature  Is  not  exponential  and  It  may  be 
known  only  numerically.  Since  the  perturbation  procedure  of  Section  V 
is  not  restricted  to  analytic  representations  of  the  streamline  curvature, 
we  used  it  to  obtain  the  correction  to  the  Gortler  number  due  the  varia¬ 
tion  specified  by  Eq.  (154)  as  an  example.  For  this  case,  the  integral 
(147a)  wa^  evaluated  numerically  while  the  integral  (147b)  was  obtained 
analytically.  Figure  7  shows  the  neutral  stability  curve  for  c  ■  0.1. 

The  Gortler  number  in  this  figure  is  based  on  the  curvature  of  the  wall 
(i.e.,  at  Hi  =  0).  Also  shown  in  Fig.  7  is  the  calculation  of  Herbert 
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using  the  Gortler  model.  We  see  that  the  decay  of  the  streamline 
curvature  stabilizes  the  flow.  The  stabilization  effect  is  pronounced 
in  the  small  wavenumber  region.  This  is  again  because  for  small  wave- 
numbers  the  disturbances  extend  farther  outside  the  boundary  layer.  As 
noted  by  Herbert^,  the  stabilizing  effect  of  the  decay  in  the  stream¬ 
line  curvature  from  the  value  at  the  wall  is  due  to  the  decrease  in  the 
driving  centrifugal  forces  associated  with  the  decreasing  curvature. 

D.  Growth  Rates  for  the  Blasius  Flow 

Curves  of  constant  growth  rate  were  obtained  for  the  Blasius  flow. 

The  results  are  plotted  in  Fig.  8,  which  gives  the  Gortler  number  G^  as 
a  function  of  the  wavenumber  a.  The  difference  between  these  curves  and 
those  obtained  by  Floryan  and  Saric  appears  only  for  o  <  0.1.  This  Is 
because  of  the  fast  decay  of  the  solution  outside  the  boundary  layer  for 
a  >  0.1.  Figures  9a-9c  show  the  eigenfunctions  u,  v,  and  w>  respectively, 
for  the  case  of  netural  stability  (i.e.,  a  =  0)  and  for  amplified  vortices 
with  o  =  2.0  at  the  wavenumber  a  =  0.5.  The  Gortler  numbers  are  1.2200 
for  o=0  and  5.4523  for  o  =  2.0.  The  persistence  of  the  disturbance 
outside  the  boundary  layer  (i.e.,  \p  >  5  where  U0  ~  0.99)  for  the  neutral 
stability  case  is  clear  from  these  figures.  Also  we  note  the  large 
gradient  of  the  solution  near  the  wall  for  o  =  2.0.  However,  no  ortho¬ 
normalization  was  needed  for  the  basic  system  while  one  orthonormaliza¬ 
tion  at  ip  =  2.328  was  needed  for  the  adjoint  system.  As  the  growth  rate 
increases,  the  number  of  orthonormalizations  Increases.  The  maximum 
number  of  orthonormalizations  was  3  for  the  results  of  Fig.  8. 


E.  The  Falkner-Skan  Profiles 


The  Falkner-Skan  profiles  correspond  to  constant  values  of  the 
pressure-gradient  parameter  Bo  defined  in  Eq.  (34).  For  flows  along 
walls  with  general  curvature  variation,  the  boundary- layer  solution  Is 
nonsimilar.  The  consideration  of  the  self-similar  solutions  In  this 
section  is  an  approximation.  That  is  we  patch  the  solution  of  the 
nonsimilar  boundary  layer  to  a  self-similar  profile  that  has  the  local 
value  of  the  pressure-gradient  parameter.  This  amounts  to  neglecting 
the  history  of  the  boundary  layer  up  to  the  position  under  consideration 
and  neglecting  the  gradient  of  Bo  at  the  local  position. 

Figure  10  shows  the  neutral  stability  curves  for  different  values 
of  the  pressure-gradient  parameter  Bo-  Favorable  pressure  gradients 
(i.e.,  Bo  >0)  are  stabilizing.  For  Bo  =  0.1,  the  Gortler  number  shows 
a  minimum  at  a  ~  0.09.  Adverse  pressure  gradients  have  a  destabilizing 
effect.  This  effect  diminishes  as  the  wavenumber  increases.  All  wave- 
numbers  smaller  than  a  certain  value  that  depends  on  Bo  appear  to  be 
unstable.  For  Bo  =  -0.05,  the  wavenumber  at  which  there  is  no  solution 
(i.e.,  no  Gortler  number)  for  neutral  stability  is  a  '  0.024  and  for 
Bo  =  -0.1  it  is  a  ~  0.048.  We  attribute  this  behavior  to  the  normal 
component  of  the  basic  flow  which  increases  as  the  pressure  gradient 
becomes  more  adverse.  By  considering  neutral  stability,  we  suppressed 
the  change  in  the  shape  of  the  vortices  in  the  downstream  direction. 

This  seems  to  be  inconsistent  with  the  growing  boundary  layer  manifested 
by  the  transverse  component  V0.  By  considering  a  parallel  stream  In  the 
basic  state  (i.e.,  using  the  Gortler  model  for  which  V0  ■  H0  =  G0  3  0 
for  all  ijJ)  and  for  Bo  =  -0.1,  we  obtained  a  curve  on  the  G^  -  a  chart 
that  is  very  similar  to  the  case  of  no  pressure  gradient  and  there  is  no 
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limiting  wavenumber  for  neutral  stability.  The  smallest  wavenumber 
considered  was  10"4.  The  value  of  G^  approaches  approximately  0.45  as 
a  -*•  0. 

Figure  11  shows  the  Gortler  number  for  neutral  stability  as  a 
function  of  Bo  in  the  range  -0.18  <  Bo  i  0.1  for  a  =  0.1.  The  separa¬ 
tion  profile  corresponds  to  Bo  =  -0.1988.  The  stabilizing  effect  of 
favorable  pressure  gradients  and  the  destabilizing  effect  of  adverse 
pressure  gradients  are  clear  from  this  figure. 

F.  Flows  over  a  Hump 

We  consider  the  curved  wall  given  by  Eq.  (38).  The  inviscid  flow 
solution  was  obtained  by  using  Davis'  method.  The  nonsimilar  boundary- 
layer  equations,  Eqs.  (31)-(33),  were  solved  by  using  a  second-order 
accurate  finite-difference  scheme  with  Aq  =  0.01  and  AS  =  0.005.  The 
asymptotic  condition  (36)  was  applied  at  n  5  8.4,  which  corresponds  to 
\p  =  11.8.  This  value  of  ^  was  taken  as  <j>m  for  the  stability  computations. 
The  numerical  values  of  the  wall  considered  are: 

H*/L*  =  G.01 ,  B  =  10.0 
The  maximum  concave  curvature  is 

k*l*  =  0.2716  at  x*/L*  =  0.845 
The  curvature  is  zero  at  x*/L*  =  0.915.  The  maximum  convex  curvature  is 
k*L*  =1.0  at  x*/L*  =  1.0. 

The  pressure-gradient  parameter  Bo  for  this  wall  is  shown  in  Fig.  12, 
while  the  curvature  variation  is  shown  in  Fig.  3. 

Figure  13  shows  the  growth  rate  a  as  a  function  of  x*/L*  for 
vortices  having  the  wavelength  X*/L*  =  0.166  and  Re  =  U*p*L*/u*  =  5.74 
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x  104.  These  values  correspond  to  neutral  stability  at  x*/L*  =  0.4.  At 
this  location  g0  3  -0.011  and  k*L*  =  0.00496.  We  see  from  Fig.  13  that 
the  vortices  grow  in  strength  in  the  downstream  direction.  The  maximum 
growth  rate  appears  at  x*/L*  =  0.75.  At  this  location  6o  3  -0.01  and 
k*L*  =  0.154.  Downstream  of  x*/L*  =  0.75  the  pressure  gradient  Is 
favorable.  The  growth  rate  drops  very  rapidly  to  zero  at  x*/L*  ~  0.83. 

At  this  location  So  =  0.283  and  k*L*  =  0.2647.  Downstream  of  x*/L*  ~ 
0.83,  the  vortices  decay  due  to  the  favorable  pressure  gradient.  Also 
shown  in  Fig.  13  are  the  growth  rates  of  the  same  vortices  calculated  by 
using  the  Blasius  profile  and  self-similar  (Fal kner-Skan)  profiles.  The 
agreement  among  these  curves  is  reasonable  in  the  region  where  the 
pressure  gradient  is  slightly  adverse.  However,  in  the  region  of  favor¬ 
able  pressure  gradient,  the  Blasius  and  Falkner-Skan  profiles  overpredict 
by  great  amount  the  growth  rates.  Thus,  using  nonsimilar  profiles  is  a 
must  in  determining  the  stability  over  such  surfaces. 

Figure  14  shows  the  results  corresponding  to  vortices  having  the 
wavelength  A*/L*  =  0.0247  when  Re  =  2.871  x  105.  These  values  were 
chosen  to  obtain  neutral  stability  at  x*/L*  =  0.4.  The  same  general 
characteristics  of  the  growth  rate  curve  of  the  previous  wavelength  are 
also  observed  in  this  case. 
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Figure  4.  Neutral  staility  For  the  Blasius  flow; 

(1)  The  Gortler  model,  (2)  The  present 
model,  (3)  Floryan-Saric  model,  (4)  The 
modified  Smith  model. 
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Figure  5.  Neutral  stability  of  the  Blasius  flow; 
effect  of  displacement  thickness. 
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7.  Neutral  stability  of  the  Blasius  flow; 
effect  of  decaying  streamline  curvature. 


Figure  9a.  Eigenfunctions  for 
the  u-component. 
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Figure  9b.  Eigenfunctions  for  the  Blasius  flow; 
the  v-component. 
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Figure  10.  Neutral  stability  for  the  Falkner-Skan  flows. 
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i 1 es  for  wavenumber 


Figure  13.  Growth  rate  of  flows  over  a  hump  for  wavelength 
A*/L*  =  0.166. 
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APPENDIX  A 


The  Navier-Stokes  Equations 
in  an  Orthogonal  Curvilinear  Coordinate  System 


Let  <j>*,  ij/*,  and  z*  denote  an  orthogonal  curvilinear  set  of 
coordinates.  This  system  is  constructed  from  the  cylindrical  surfaces 


4>*(x*,y*)  =  C^ 

(Al) 

4>*(x*,y*)  = 

(A2) 

and  the  planes 

z*  =  C2 

(A3) 

where  C^,  C^,  Cz  are  constants. 

We  consider  a  steady  incompressible 

flow  with  constant  properties.  The  equation  of  motion  in  the  coordinate 
system  ( A1 ) - ( A3 )  are19 


4)*-Momentum  Equation 


y_3y.+  v_3u.+  3u  /_v _u ah(j>x  _  1_11E 

%  %  9*  32  "  hfy  d*  hfy  ^  V  3 * 

+  h  «  32u  +  h  92w  .  A  i_  il^\  .  jfc  3  ,u  n 

h4>n^  dz7  ^  9$  ^h^  h^  ^h^  ’ 


i  ■>„  ,,  3h.  3h, 

hA  *  }  ** 


(A4) 
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4*-Momentum  Equation 


L  f\L 

y_9v  v_9v  av,  ,_v _ _jJ>  u  _  4,  _  1_  1  Ifi 

\  84>  \  94  92  \\  9$  '  ^  ^  P  3* 

k2  U 

,  V  3  r  4  9  /V\j_u  9  /  ^Mj.9  fo  4  9v  a  ,j  U  4\ 

vv  *  cv  *  V  ♦  *  ^  ^  ^ 3* 


4 


4 


<  •  i 

4 


+  h  s2w  +  h  h  ^X  +  —  d—)  +  r^-  —  (— )]  • 

"4  343z  n<p\  dz7  Lh^  34  34  34 


1  ^  9h,  9h, 

p  /L.  9y  +  .  v. _ 4y  ._£ 

'h,  34  h.h,  34  '  34 

4  4  4 


(A5) 


z*-Momentum  Equation 


u  3w  v_  3w  ,  3w  _  X  9£.  +  v  3  /  l  3u^  ,  _jj>  3W\ 

34  34  W  3t  "  "  p  3z  34  3z  h^  3 <fiJ 

h  ■> 

,  3  /  d>  3w  ,  ,  3v  \  ,  9u  l.  3  w 

34  34  h4  3z  2h4h4  9z^ 


(A6) 


Continuity  Equation 

If  <y> *  k  <y> +  y*  If  - 0  <A7> 

All  the  variables  in  the  above  equations  are  dimensional;  we  dropped  the 
asterisks  for  simpler  notation.  Where  u*,  v*  and  w*  are  the  velocity 
components  in  the  4*,  4*,  and  z*  directions  respectively;  p*,  p*  and  v* 
are  the  static  pressure,  density,  and  kinematic  viscosity  respectively; 
h*  and  h*  are  the  metric  coefficients  for  the  curves  yiven  by  Eq.  (A1 ) 
and  (A2). 

The  Navier-Stokes  equations  when  4*  and  4*  are  the  potential  and 
stream  functions  of  a  two-dimensional  incompressible  irrotational  flow 
become 
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<t*-Momentum  Equation 


hu 


au 

34 


+  hv 


3u 

34 


h2 


w  -  v*  +  UV 


3  Z 


Ml 

3$ 


-|r  =  -  -  Hr  +  v[72u 
34  p  34 


32h\  ,  2  /  3h  3v  3h  3v \ -i 
3ip';  h  v34  34>  "  34>  3i pn 


^♦-Momentum  Equation 


3  v 


hu  It  +  hv  U  +  h*W  32 


-  li2 


3h 

34 


32h\  ,  2_  / 3h  3u  3h  3u \ ~i 
dip*'  h  ^4  dip  “  dip  d<f>  ^ 


2*-Homentum  Equation 


.  3w  , 

hu  3*  + 


hv 


3w 

34 


+  h!«  H  =  - 


hi 

p 


If  +  vV2w 


Continuity  Equation 


3hu  ,  3hv  ,  l,2  3w  _  n 
W  W  h  32  "  ° 


where 


h  = 


V2 
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Appendix  B 

Schwarz-Christoffel  Transformations  for  Curved  Surfaces 


Davis^4  presented  a  numerical  method  for  integrating  Schwarz- 
Christoffel  transf ormations  for  curved  surfaces.  Sufficient  information 
about  the  method  and  several  examples  are  given  by  Davis. 

Here,  we  specialize  his  method  for  the  wall  given  by  Eq.  (38), 
that  is 

y  =  H  sechB(1  •  •;)  for  |x|  <  *  (Bl) 

We  assume  that  the  wall  is  specified  by  Eq.  (Bl)  in  the  interval 
0  <  x  <  2  and  it  coincides  with  the  x-axis  outside  that  interval.  For 
the  case  treated  in  Section  II,  H  =  0.01  and  B  =  10,  hence 

y(0)/H  =  y(2)/H  =  2e'10  «  1 

Therefore  the  above  assumption  is  justified. 

We  divide  the  wall  in  the  interval  0  <  x  <  2  into  a  finite  number 
of  elements  N.  Let 


z„  =  xm  +  iy„ 
m  m 


(B2) 


denote  the  left-end  point  of  the  inth  element, 


(B3) 


denote  the  angle  between  the  tangent  line  and  the  x-axis  (positive 


clockwise)  at  zm,  and 


ew  -  -  tan",(x 


i  /^m+1  •'m 


"  Yrr 


) 


(B4) 


m+1  m 

denote  the  angle  between  the  straight  line  connecting  the  mth  and 

(m  +  1 ) th  points  on  the  wall  and  the  x-axis  (positive  clockwise).  As 
14 

shewn  by  Davis  ,  a  transformation  which  maps  the  upper  half  of  the 
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5-plane  onto  the  region  above  the  curved  wall  in  the  z-plane  is 

dz  m  n  C2rn/7T  C3m/n 
d5  =  M  g2m  93m 


(B5) 


where 


%  * 


c3m  .  .  em)/(bm+1  - 


(B6) 


m+]  ~ 
bm+l  ”  bm 


(b, 


m+1 


bm>C3m 


(B6) 


9 2m  =  Am/Am+1 


(B7) 


^3m  ^n/^m+1 

U  -  b  )  b 

*m  ‘  -  bm)  m  e  m 


B  =  (5  -  bm) 
m  ^  m' 


(c2  -  ti)  u  +  \)2n 


(B8) 

(B9) 

(BIO) 


Here  the  bm  are  images  of  the  zm  in  the  5-plane;  they  are  real  numbers. 

The  unknowns  in  Eq.  (B5)  are  M  and  the  bm.  They  are  determined 
iteratively  as  follows: 

(a)  We  guess  values  for  the  bffl  and  M.  In  the  present  applica¬ 
tion  M  is  purely  real;  however,  this  need  not  be  specified 
since  it  will  come  from  the  solution  procedure  automatically. 

The  bm  are  picked  such  that  b_  ,  <  b  <  b, , , . 
m  m- 1  m  m+1 

(b)  By  using  complex  arithmetic,  we  integrate  Eq.  (B5).  We  set 
bi  =  0  (i.e.,  the  point  z  =  0  is  transformed  to  5  *  0).  Now 
the  obtained  Zgm  will  not  coincide  with  the  actual  zm  given  by 
Eq.  (B2).  A  correction  to  M  Is  obtained  first  by  requiring 
that  the  point  z  =  (1.0,  H)  be  transformed  to  the  point 
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C  =  (1.0,  0).  Updated  values  for  the  bm  are  obtained  from 


bcm  '  bcm-l  _  ,  lzm  zm-l 

b  -  b  “  K 
gm  gm- 1 


zgm  ‘  zgm-l 


(Bll) 


where  the  subscript  c  denotes  corrected  or  updated  values  and 
the  subscript  g  denotes  guessed  or  old  values.  The  scale  k  Is 
related  to  M  by 


c.  A  new  iteration  is  performed  by  using  the  new  values  of 
M  and  the  bm  (step  a).  Step  (b)  is  repeated  until  conver¬ 
gence  is  achieved. 

The  complex  velocity  in  the  z- plane  is  obtained  as  follows.  Let 

af  *  (biz) 


In  the  c-plane  the  flow  is  a  uniform  stream  parallel  to  the  £-axis. 
Hence 


(B13) 


where  w  is  the  complex  potential. 

(u  .  iv)  =  dw  =  dw 
'  v  z  dz  d?  dz/d£ 


In  the  z-plane 


(u  "  iv)z  =  Rf  (B14) 

N 

Since  l  6^  =  0  then  f  ->  1  as  c  -*■  «.  The  uniform  flow  in  the  c-plane 
m=l 

is  mapped  onto  a  uniform  flow  in  the  z-plane.  Hence 


(B15) 
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Substituting  Eq.  (B15)  into  Eq.  (B14)  we  obtain 


u  -  iv  =  f  (B16) 

In  the  actual  calculations  U  =  1 .  The  metric  coefficient  h  Is  thus 
given  by 


h  =  |f  | 


(B17) 


f 
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Appendix  C 
The  Matrix  AA 


The  matrix  AA  is  a  6  x  6  square  matrix.  The  non-zero  elements  of 
this  matrix  are 

a2i  =  (h2  -  1  )ct2  -  V0hK  -  ^  ~  (°2  +  ^  So  -  Soa  -  $o  ^jp) 

^  HrK 


(Cl) 

323  =  |tt  (So Vo  +  2ahK^  -  Bo^^)  (C2) 

324  -  ^  (C3) 

326  =  (C4) 

3  3  3  =  hK^  ( C5 ) 

a34  =  -  (h  -  l)a  (C6) 

a54  =  (h2  -  1  )a2  -  £  (C7) 

d<j> 

356  =  -  (h  -  l)o  (C8) 

36 3  =  \  SohK^  +  ahK^  +  \  SJ  -  2Uo(h  -  1 )G2  (C9) 


a&3  =  -(h2  -  1  )a2  2V„hKj>  +  (hty  *  +  h2K2  ♦  £*  |  ^ 


(CIO) 

364  =  v0(h  -  1 )a 

(Cll) 

36  5  =  ■  ( h  -  1  )ct 

(C12) 

366  =  -hK^ 

(C13) 

The  prime  denotes  differentiation  with  respect  to  ip. 
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Appendix  D 

The  Matrices  P  and  P* 

These  matrices  are  given  in  the  next  two  pages. 
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